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Abstract. - We show how inelastic collisions induce a power law with exponent —3 in 
the decay of the angular velocity distribution of anisotropic particles with sufficiently small 
moment of inertia. We investigate this question within the Boltzmann kinetic theory for an 
elongated granular particle immersed in a bath. The power law persists so long as the collisions 
are inelastic for a large range of angular velocities provided the mass ratio of the anisotropic 
particle and the bath particles remains small. Suggestions for observing this peculiar feature 
are made. 



Granular gases are systems in which macroscopic particles lose a fraction of their kinetic 
energy at each collision [1,2]. In the absence of energy supply, starting from an homogeneous 
state, the granular fluid cools down [3] and after a finite time one observes a spontaneous 
symmetry breaking with the formation of clusters as well as of shear waves and convection. 
Conversely, if an external energy supply is continuously brought to the particles, the sys- 
tem may reach a non-equilibrium steady state (NESS), whose properties differ significantly 
from those of thermal equilibrium (breakdown of the equipartition [4,5], non-gaussian statis- 
tics, modified hydrodynamics [6],...). All those characteristics are intimately related to the 
dissipative nature of collisions. 

The studies of granular gases have focused mainly on spherical particles [7] . Recent work [8] 
on nearly smooth spherical grains stressed the relevance of rotational degrees of freedom for 
the hydrodynamic behavior. The rotational energy is ubiquitous for anisotropic particles, and 
one can reasonably expect that the anisotropy (generally present in real granular systems) 
might introduce additional effects. Only few studies have addressed up to now this question 
(free cooling state in a three-dimensional system of needles [9] , breakdown of the equipartition 
between different degrees of freedom [10,11]. Experimentally, the dynamics of shaken granular 
dimer gas has been investigated [12]. 

In this letter, we show that the stationary angular velocity distribution of a thin granular 
needle of mass M acted upon by inelastic collisions in a bath of point masses m shows a power- 
law decay in a large range of angular velocities provided M/m <^ 1. We relate this feature to 
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the properties of a rigorous solution of the Bohzmann equation in the Mjra hmit.The 
region Mjm ^ 1 studied here is opposite to the Brownian motion case Mj m ^ 1 . The 
latter has been thoroughly studied for granular fluids at least for spherical particles [13,14]. 
We discovered that the dynamics of a low mass anisotropic particle is qualitatively different. 
It turns out that when the ratio M/m is finite but small this unphysical (zero moment of 
inertia) solution continues to manifest itself through a power-law decay observed in a large 
range of angular velocities whereas a gaussian-like decay is restored only for the high-energy 
tail. When the restitution coefficient is set to 1, this peculiar behavior disappears. Clearly, 
in real systems collisions with the ends of the elongated particle should be also considered. 
The analysis performed in Ref. [11] has shown that when the length of the particle is large 
compared to its width and also to the size of the bath particles the coUisional contributions 
from the extremities are negligible and predictions of the simplified model considered here 
should apply. We also present suggestions for an experiment which would permit to observe 
the power-law decay of the angular distribution function. 

Quite recently a class of stationary states with power-law high-energy tails (implying an 
infinite energy) has been found for homogeneous granular systems with the suggestion of 
implementing them by energy injection at large scales only [15,16]. It should be stressed 
however, that the underlying mechanism is completely different from the one found here. We 
deal with an impurity (the needle) immersed in the granular fluid (which is in a NESS). The 
power law appears owing to the coUisional rescaling of the distribution of the bath with the 
weight depending on the position of the impact point. The anisotropy of the impurity plays 
in this rescaling an important role. The possibility of a power law for an impurity has been 
already mentioned in the review on the Maxwell model [17]. However, to get this effect the 
same power law had to be assumed for the granular bath. Our results reveal an essentially 
new mechanism for a physically relevant case of a free motion between collisions. 

We consider a two-dimensional system consisting of an infinitely thin homogeneous needle 
of length L and moment of inertia / = ML^/12, with the fixed center of mass. The needle 
undergoes inelastic collisions with the bath particles. Between collisions it rotates freely 
around the axis passing through its center and perpendicular to the plane of motion. The 
only degree of freedom of the needle is its orientation specified by a unit vector u that points 
along its axis. The rate of change of the orientation u = uj\i±_ involves the angular velocity 
we] — oo, +oo[ and a unit vector u_l perpendicular to u. 

At a binary collision the needle is hit by a point mass of the bath at some point Au, |A| < 
L/2 where A is the algebraic abscissa along the needle (see Fig^). The relative velocity V at 
the point of impact equals 



where v denotes the velocity of the bath particle. 

The instantaneous binary collisions conserve the total angular momentum. Labelling the 
post-coUisional quantities with a star and using for any vector w the convenient notation 
w±_ = w • u^, w\\ = w • u we thus write 



V = v — Au = v — Awuj^ 



(1) 



loj* 4- Xmv*^ = luj + Xmv± 



(2) 



The relative velocity will be assumed to change according to the collision law 




VI = -~aV± 



(3) 
(4) 



involving the normal restitution coefficient < a < 1. 
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By combining Eqs. C3)-I0J, we find the collisional change of the particle velocity 

v±=vj_ — — — . (5) 

1 + mA^ 



and the corresponding change of angular velocity 

^ = ^ + ^ r X9 • 



(6) 



The inverse transformation is obtained by replacing a by a~^. 

The stationary Boltzmann equation for the angular velocity distribution ^"(0;) of the nee- 
dle expresses the invariance of the distribution under collisional processes. The gain term 
corresponding to the post-coUisional angular velocity lo must be exactly compensated by the 
loss term involving the angular velocity lo as the pre-coUisional one. The equation reads 

dX J d^\vi_ - Xlu\ l^^^f-J-p L - F(a;)$B(v)j = 0, (7) 

where the pre-coUisionnal velocities v** and uj** follow from Eqs.©-® by switching the 
roles of initial and final velocities and replacing a by a~^. $b(v) is a velocity distribution 
describing the steady state of the bath. 

The integration over the velocity component v^^ in Eq.(|7I) can be readily performed as this 
variable does not show in the arguments of F. Putting then vj_ ~ {XLu)y one finds the integral 
equation 



'L/2 



dXX' J dy\y\F (^co (^1 + y^-^ 



(l+a)mX'^\\ , / / (amX^-I) 
^ ^ (l>B [ Xlj il + y' 



+ mA^ / / V V / + mX 



= F{lu) f ' dXX' f dy\y-l\4>B{XLoy) (8) 

J-L/2 J 

with (I)b{v) = / dv||$B(|v|) 

We consider first the limit where both the mass of the needle and the coefficient of restitu- 
tion a tend to zero. In this case the velocity of the bath particles is not modified by collisions 
whereas the angular velocity of the asymptotically massless needle is reset after each collision 
acquiring instantaneously the value u±/A, i.e. the ratio of the bath particle approach velocity 
to the abscissa of the impact point. The Boltzmann equation ((HI becomes 

dAA2 / dy\y\F {uj (1 + y)) <^b {Xoj) = Fiuj) [ dXX^ [ dy\y ~ l\M>^^v) (9) 

L/2 J J J 

The exact solution F{uj) of Eq.© is obtained by scaling the bath distribution (j>B{v) with the 
weight depending on the position of the impact point 

F{cj)= J ^^^dXi^ — j ^BiXiv) (10) 

It should be stressed that this result does not require specific assumptions about the distri- 
bution of the bath particles. Indeed, Ea.l(Tn|l imphes the power law decay F{lu) oc lo^^ for 
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Fig. 1 - The needle and a bath particle in the plane: u is a unit vector along the axis of the needle 
forming with an orthonormal system. The element of the needle at Au moves with linear velocity 
Aa;u^ . For a point of impact Au the pre-coUisional relative velocity equals V = v — \ljuj_ 

large values of w provided the second moment of the distribution 0s (w) exists. The specific 
shape of (j)B{v) is irrelevant. We also note that for the rescaled distribution ((117)) the second 
and higher moments are not defined. It is thus natural to think that the stochastic variable 
u) is subject to Levy flights. 

When the bath is at thermal equilibrium described by the Maxwell distribution (/>b(v) = 
\/m/ 2'kT exp ( — mw^ / 2T) , the stationary distribution takes the form 

4 / mL^u^\ 

The power-law is reached for uj > \/8T/ (mL^). Fig. [3 shows the logarithm of the distribution 
function versus the angular velocity. The inset in Fig|2 displays the log- log plot illustrating 
the crossover to the power-law behavior. 

Eas. 1)10111(1 correspond to the maximal inelasticity a = 0. When 7\/ = and < a < 1, 
the solution of the Boltzmann equation behaves asymptotically as 
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Fig. 2 - Log-linear plot F{ijj) of the "zero mass" needle. The inset displays the log-log plot showing 
the power-law decay for sufficiently large lo values for T/{mL^) — 1/2. 

Direct numerical integration of Eq.JSJ (with 7 = 0), along the same lines as in [18], is shown 
in Figl^lfor various values of a. The inset of Fig |2| displays H{lj) = F(cj) versus 
LO, indicating that the asymptotic behavior is rapidly reached. It should be noticed that for 
elastic collisions (a = 1) all these effects disappear as then 

(13) 

and the limit / ^ does not yield a probability distribution. 

We turn now to the physically relevant case of a needle with a small but finite mass 
< Af ^ m. An asymptotic analysis of the distribution F{u!) can be then performed on the 
basis of Eq.®. Assuming that at large angular velocities the distribution function is gaussian, 
F{'w) ^ cxp(— /cj^/ (2T) (an assertion well supported by numerical results [19]), one can show 
that r = (1 + a)T/2, irrespective of the mass ratio of the needle and the bath particles (it 
is an exact asymptotics in the Brownian motion limit M/m — > oo). By means of an accurate 
numerical solution of the Boltzmann equation (details will be given elsewhere [19]), we could 
essentially confirm this asymptotic behavior inferred first by analytic arguments. In addition, 
in the case under consideration (0 < M /m <C 1), we could identify a sub-leading multiplicative 
term (decreasing less rapidlly than a gaussian) which depends on the mass ratio as well as on 
the coefficient of restitution. This combined behavior occurs in the asymptotic regime. 

However, our really important conclusion is that before the very high energy asymptotics 
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Fig. 3 - Log-Log plot of F{uj) of the "zero mass" needle with different values of the coefficient of 
restitution a — 0, 0.2, 0.4, 0.4, 0.8 from top to bottom. The inset shows the oj-dependence of the 
rescaled function H{lj) = F(a;) o;^ for T/{mL^) = 1/2. 



is attained the power-law ~ uj^^ found for a massless needle remains present with a crossover 
occuring only at uJc ~ ^J2T/I. This situation is illustrated in Fig ^ where the log- linear plot 
of the angular distribution functions is displayed for mass ratios M/m = 0,0.005,0.01, and 
the dashed curve represents the zero- mass limit. For lj < ujc, one observes the power-law 
decay (see the log-log plot in the inset) whereas the gaussian behavior is recovered beyond 
(jOc- Contrary to the zero mass case, the existence of a gaussian decay at very large values 
of the velocity leads to a finite granular temperature for a needle with a finite mass, and 
more generally all moments of the distribution functions are now well defined. In fact, one 
can show that the granular temperature defined as the second moment of the distribution 
times the moment of inertia goes to zero with the mass of the needle if a < 1, whereas the 
temperature remains independent of the mass if a = 1. This result is also a consequence of 
the existence of a well defined solution of the Boltzmann equation for a zero-mass needle when 
a < 1. In conclusion, the behavior for M/m ^ 1 is qualitatively different from that observed 
in the Brownian motion limit involving an new coUisional scaling mechanism producing the 
w"^ power law. 

Since the power-law regime corresponds to low and moderate angular velocities, it should 
be accessible to an experimental observation. Recent experiments on intensely vibrated gran- 
ular systems using high speed photography [20,21], image analysis techniques and particle 
tracking [22,23] show that many quantities (granular temperature, velocity profiles,...) can 
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Fig. 4 " Log-Log plot of F{uj) for a needle with a mass M/m = 0.005, 0.01. (dashed curve corresponds 
to M = 0.0). The inset shows the crossover between the power law and the gaussian-like asymptotics 

be precisley measured allowing the study of the equipartition or fluctuation theorem. Several 
orders of magnitude of velocities can be monitored [12,24] within the region where the power- 
law decay is expected. More specifically, by immersing in a two-dimensional granular gas a 
thin rigid needle which rotates around a perpendicular axis, it should be possible to obtain 
the angular velocity distribution in the relevant region of lo. By using different needle lengths 
or different sizes of the bath particles one should observe successive crossover velocities. We 
expect that the collisional rescaling found here will induce the power law also in the case of 
anisotropic particles of a general convex shape in the region accessible to experiment. 

J. P. acknowledges financial support by the CNRS, France, and the hospitality at the Lah- 
oratoire de Physique Theorique de la Matiere Condensee, UPMC (Paris) where this research 
has been carried out. 
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